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Abstract 

Using the dielectric theory for a weakly coupled plasma we investigate 
the stopping power of the ion in a temperature anisotropic magnetized 
electron plasma. The analysis is based on the assumption that the energy 
variation of the ion is much less than its kinetic energy. The obtained 
general expression for stopping power is analyzed for weak and strong 
magnetic fields (i.e., for the electron cyclotron frequency less than and 
greater than the plasma frequency), and for low and high ion velocities. It 
is found that the friction coefficient contains, in addition to the usual ve- 
locity independent friction coefficient, an anomulous term which diverges 
logarithmically as the projectile velocity approaches zero. The physical 
origin of this anomulous term is the coupling between the cyclotron mo- 
tion of the electrons and the long-wave length, low-frequency fluctuations 
produced by the projectile ion. 
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I INTRODUCTION 



Energy loss of the ions in a plasma has been a topic of great interest due to its considerable 
importance for the study of basic interactions of the charged particles in real media. Recent 
applications are electron cooling of heavy ion beams |], |3|] and energy transfer for inertial 
confinement fusion (ICF) (see || for an overview). Electron cooling is realized by mixing 
the ion beam periodically with a cold electron beam of the same average velocity. The 
interaction length is normally about a few meters and the electron beam is guided by a 
magnetic field parallel to its direction of motion. The cooling of the ion beam may then 
be viewed as an energy loss in the common rest frame of both beams. Similar questions 
arise in heavy-ion-induced ICF. There a frozen hydrogen pellet is heated and compressed by 
stopping of ion beams in the surrounding converter. In this case the electrons of the solid 
state converter are acting like a plasma and absorb the incoming energy. 

In the electron cooling process the velocity distribution of the electron beam is highly 
anisotropic because of the acceleration from the cathode to the cooling section. It can be 
described by a Maxwell distribution with two different temperatures, a longitudinal Tji and 
a transversal T± \q t|. Furthermore, an external, longitudinal magnetic field is needed to 
guide the electrons from the cathode to and through the electron cooler and to stabilize the 
anisotropic velocity distribution by suppressing the transverse-longitudinal relaxation. 

In the present paper we are interested in the influences of the magnetic field and the 
temperature anisotropy on the ion beam stopping power. 

Since the early 1960's several theoretical calculations of the stopping power in a mag- 
netized plasma have been presented |], [7| 0, [10], [Ll], 0, 0- Stopping of a fast 
test particle moving with velocity V much higher than the electron thermal velocity v t h 
was studied in Refs. [|, |6], Energy loss of a charged particle moving with arbitrary 
velocity was studied in Ref. 0. The expression obtained there for the Coulomb logarithm, 
A = \n(\[)/pj_) (where A^> is the Debye length and p± is the impact parameter for scat- 
tering for an angle $ = vr/2), corresponds to the classical description of collisions. In the 
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quantum-mechanical case, the Coulomb logarithm is A = ln(Ac/As), where As is the de 



Broglie wavelength of plasma electrons [15 



In Ref. ||10|| , the expressions were derived describing the stopping power of a charged 
particle in Maxwellian plasma placed in a classically strong (but not quantizing) magnetic 
field (Ab C a c C A^, where a c is the electron Larmor radius), under the conditions when 
scattering must be described quantum mechanically. Calculations were carried out for slow 
test particles whose velocities satisfy the conditions (m/mj) 1 / 3 f t h < V <C v t h, where m ; is 
the mass of the plasma ions and m is the electron mass. 

In the recent paper (jTTJ the stopping power in the magnetized plasma has been inves- 



tigated for high-velocity light particles taking into account the Larmor rotation of a test 
projectile in a magnetic field. It has been shown that the stopping power can exhibit an 
oscillatory dependence on the magnetic field and that it is much greater than in the case 
without magnetic field. 

More attention has been paid on the stopping power in a strongly magnetized plasma 
for ions which move along the magnetic field [11, [12|, 13|. Both uncorrelated Jll], [13] and 



correlated |12j situations have been discussed. 

These investigations have concentrated on the stopping power in temperature isotropic 
plasma. Extensions to nonlinear effects of ion stopping and temperature anisotropy have 
been done recently by particle-in-cell (PIC) computer simulation where the case Tj| <C 
T± has been investigated which is interesting for electron cooling process. Here, in the 
framework of dielectric theory, we will focus on the stopping power at arbitrary temperature 
anisotropy T±/T\\. 

The paper is organized as follows. We start in Sec. II, with solving the linearized Vlasov- 
Poisson equations by means of Fourier transformation. This provides the general form of the 
linearized potential generated in a temperature anisotropic magnetized Maxwellian plasma 
by a projectile ion from which the stopping power is deduced. 

In the next Sec. Ill, is dedicated to apply our results to nonmagnetized plasma. Cal- 
culations are carried out for small projectile velocities at arbitrary temperature anisotropy 



and arbitrary direction of ion motion with respect to the anisotropy axis. 

Then we turn to the effect of a weak magnetic field on the stopping power in Sec. IV, 
while we concentrate on the influence of a strong magnetic field in Sec. V. In contrast with 



the papers |TT], |T^] we consider an ion motion in arbitrary direction. 

As the last issue we investigate in Sec. VI the stopping power for small projectile veloc- 
ities at arbitrary magnetic field and temperature anisotropy. The friction coefficient there 
contains an anomalous term which increases logarithmically when the projectile velocity 
approaches to zero. 

The achieved results are finally summarized and discussed in Sec. VII. 



II DIELECTRIC THEORY 

For the temperature anisotropic plasma with two different temperatures Tjj, T± of the 
electrons we define an average temperature T = |Tj| + §Tj_. Within the dielectric theory the 
electron plasma is described as a continuous, polarizable fluid (medium), which is represented 
by the phase-space density of the electrons /(r, v, t). Here, only a mean-field interaction 
between the electrons is considered and hard collisions are neglected and the evolution of 
the distribution function /(r, v, t) is determined by the Vlasov-Poisson equation is valid for 
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weakly coupled plasmas where the number of electrons in the Debye sphere = 47rn A D 3> 
1 is very large. Here no is the electron density, Ad = (k B T / A-Kn^e 2 ) 1 / 2 is an averaged Debye 
length. 

In the following, we consider a nonrelativistic projectile ion with charge Ze and with a 
velocity V that moves in a magnetized temperature anisotropic plasma at an angle •d with 
respect to the magnetic field B . The axis defined by B also coincides with the degree 
of freedom with temperature Tji. We assume that the energy variation of the ion is much 
smaller than its kinetic energy. The strength of the coupling betweeen an ion moving with 
velocity V and the electron plasma is given by the coupling parameter 
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Here t> t h = (k^T/m) 1 ^ 2 is the average thermal velocity of an electron. The derivation of 
Eq. (1) is discussed in detail in Ref. [ 16| . The parameter Z characterizes the ion-target 



coupling, where Z <C 1 corresponds to weak, almost linear coupling and Z > 1 to strong, 
nonlinear coupling. 

For a sufficiently small perturbation (2 < 1) the linearized Vlasov equation of the 
plasma may be written as 

dfx .dfx dfx e df 

+ v "a v X b — = j- 2 

at or av m or av 

where f — fo + fx and the self-consistent electrostatic potential cp is determined by the 

Poisson equation 

v 2 (p = -AixZe5(r - Vt) + 47re J dvfx{r, v, t). (3) 

The b is the unit vector parallel to B , — e and uo c = eB /mc are the charge and Larmor 
frequency of plasma electrons respectively, fo is the unperturbed distribution function of 
plasma electrons, which in the case of temperature anisotropic, homogeneous electron plasma 
is given by two Maxwellians for the longitudinal and transversal degrees of freedom 

Mv " v±) = w^fci exp exp 1 1 ■ (4) 

where (t|) = = k B T\\/m, (v 2 ± ) = 2v 2 h± = 2k B Tj_/m. 

By solving Eqs. (2) and (3) in space-time Fourier components, we obtain the electrostatic 
potential 

Ze r „ exp \ik(r — Vt)] 
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which provides the dynamical response of the temperature anisotropic plasma to the motion 
of the projectile ion in the presence of the external magnetic field. Here e(k, uo) is the 
dielectric function of a temperature anisotropic, magnetized plasma which is given by 



e(k,u) = 1 



with 



1 + 



PA 2 



[G(s)+iF(s)} 



1 + is\p2 / dt exp \istV2 - X(t) 



exp 



(6) 



+ - sin a{± — T) at sm — 

uj c J \ kv t h\\ , 



istV2-X(t) 



OJ, 



cos 



kv 



thi i 



(7) 

v th± /u c and 



X(t) = t 2 cos 2 a + k 2 al ± sin 2 a 

where Any = ^th||/^ P , u p is the plasma frequency, s = tu/kv^, r = T±/T\\, a c 
a is the angle between the wave vector k and the magnetic field. 

As shown in Appendix A, Eqs. (6) and (7) are identical with the Bessel function rep- 
resentation of e(k, u) derived e.g. by Ichimaru |17| . Eqs. (6) and (7) are, however, more 
convenient when studying the weak and strong magnetic field limits in Sees. IV and V. 

The stopping power S of an ion is defined as the energy loss of the ion in a unit length 
due to interactions with the plasma electrons. From Eq. (5) it is straightforward to calculate 
the electric field E = — \/ip, and the stopping force acting on the ion. Then, the stopping 
power of the projectile ion becomes 



S = — 



dE 
df 



or 



(8) 



r=Vt 



2Z 2 e 2 X 2 r 



D|| 



1 7 

k 3 dk J dfx j d(p- 





cos9F(s) 



•/ J J ' [h*\l {] +G(s)] 2 + F*(s) ' 

where fi = cos a was the angle between k and Bo, is the angle between k and V, s = 
k-V//cf t h|| = (y/ v th\\) cos0, cos 6 = fi cos $ — a/1 — H 2 sin $ cos (p, and is the angle between 
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V and B . In Eq. (8) we introduced a cutoff parameter A; max = l/r min (where r min is the 
effective minimum impact parameter) in order to avoid the logarithmic divergence at large 
k. This divergence corresponds to the incapability of the linearized Vlasov theory to treat 
close encounters between the projectile ion and the plasma electrons properly. For r min 
we thus use the effective minimum impact parameter of classical binary Coulomb collisions 
r min = Ze 2 /mv 2 for relative velocities v r ~ (V 2 + v 2 h ) 1 ^ 2 , which is often called the "distance 
of closest approach." Hence 

k - J- - m{ y 2+ ^ ( 9) 

A two temperature description of an electron plasma is valid only when the ion beam- 
plasma interaction time is less than the relaxation time between the two temperatures, T\\ and 
Tj_. For an estimate we will briefly consider the field-free case, because the external magnetic 
field suppresses the relaxation between the transversal and longitudinal temperatures during 
the time of flight of the ion beam through plasma. 

The problem of a temperature relaxation in a temperature anisotropic plasma with and 
without of an external magnetic field was considered by Ichimaru [fL7[ |. Within the dominant- 
term approximation the relaxation time Ar re i for the plasma without magnetic field is given 
by 

1 8 ^ ^ lnA c , (10) 



Ar rcl 15 Vm (fc B T eff ) 3 / 2 
where lnA c = ln(A^) is the Coulomb logarithm and the effective electron temperature T e g 
is defined through 
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3/2 



15 

y 



/i 2 (l - fi 2 )dfi 



-eff o 

5VS (l + 2r) 3 / 2 



12T 3/2 



l) s 



[» 2 T\\ + (1 
r + 2 



/i2)T ± ]3/2 



(11) 



T - II 



P0{T) 
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In l±iO=z T < i 

Po(r)H ^ . (12) 

arctan Vt — 1 , r > 1 

The relaxation time calculated from Eq. (11) are of the order of 10~ 6 s, 0.5 x 10~ 5 s and 

lCT 3 s for averaged temperatures T = 10~ 2 eV, T = O.leV and T = leV, respectively, for 

anisotropies r ~ 0.01 — 100. The interaction time (for instance, for ICF or for electron 

cooling) is about 10 -7 — 10~ 8 s. Therefore, ion beam-plasma interaction time can be very 

small compared to the plasma relaxation time. 



Ill STOPPING POWER IN PLASMA WITHOUT MAGNETIC FIELD 

Let us analyse expression (8) in the case when a projectile ion moves in a temperature 
anisotropic plasma without magnetic field. The plasma dielectric function from Eqs. (6) 
and (7) now takes the form 

^) = 1 + F^>Q) ' < 13) 

Here A = (/i 2 + r(l — /i 2 )) 1 / 2 and W(s) = go(s) + ifo(s) is the plasma dispersion function 

m 



g (s) = 1 - sV2 Di {^jT^j ! /o( s ) = \f^ s ex P 



s 



(14) 



where 



Di(s) = exp(-s 2 ) J dtexp(t 2 ) (15) 
o 

is the Dawson integral |18| which has for large arguments s the asymptotic Di(s) ~ l/2s + 
l/4s 3 . 

Substituting Eq. (13) into Eq. (8) and performing the ^-integration we obtain 
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Z 2 e 2 



Sn 



^ A D|| 



(/// / d(f^^- Qq 



v cos 6 
v t h\\ A 



(16) 



where £|| — fc max AD|| and 



Qo{x,0 = f (x) In 
9o(x) 



arctan ■ 



e+9o(x) 



(17) 



arctan 



fo(x) f (x) 

In the case of temperature isotropic plasma (T± = Tji = T, and r 
(16) coincides with the result of e.g. Ref. |19 



1) A = 1 and Eq. 



2 2 2 V/?,th 



(18) 



where v th = v th \\ = v th± , \ D = v th /u p , and ^ = k max \ D . 

When a projectile ion moves slowly through a plasma, the electrons have much time to 
experience the ion attractive potential. They are accelerated towards the ion, but when 
they reach its trajectory the ion has already moved forward a little bit. Hence, we expect 
an increased density of electrons at some place in the trail of the ion. This negative charge 
density pulls back the positive ion and gives rise to the stopping. This drag force is of 
particular interest for the electron cooling process. In the limit of small velocities 5* ~ R ■ V. 
This looks like the friction law of a viscous fluid, and accordingly R is called the friction 
coefficient. However, in the case of an ideal plasma it should be noted that this law does 
not depend on the plasma viscosity and is not a consequence of electron-electron collisions 
which are neglected in the Vlasov equation. 

The Taylor expansion of Eq. (16) for small V(V v t h) yields the friction law 



(e 2 A 



where £ = k max X D 



So — 

3V27T l>th 

(1 + V*/vl) /Z ~ l/Z 



V [/i(r)+/ 2 (r) sin 2 tf 



(19) 
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h{r) 



3 (2r + l\^ } d M^AM) 



h(r) 



and the function ip is 



3 (2r+l^ } S^M^/i)) 



2^(0 



(20) 



(21) 



V>(0 = k(l + £ a )- 



? 2 



i+e 



(22) 



In the case of temperature isotropic plasma (r = 1) we have I\ — 1 and J2 = 0. Then 
the Eq. (19) becomes the usual friction law in an isotropic plasma fll9| . For the strongly 
temperature anisotropic case, when r <C 1 (T± <C Tj|) we have £|| ~ \fZjZ and 



/1 



v/3 



Lz 2 (l + 6f) + ln(l+rf 



>/3 



£j + 21n(l + fJ)+3Li 2 (l + £j 



(23) 



(24) 



12^(0 

Here the functions 1\ and 7 2 do not depend on r, and Li 2 (x) is the dilogarithm function 
2TJ. Note that Z < 1 and therefore £ > 1, £y > 1 in Eqs. (23) and (24). The Coulomb 



logarithms in Eqs. (23) and (24) are then the leading terms and 



r A 1 ^ t v^3 1 

ii ~ — In — <s; lo — 

1 6 Z 2 8Z 2 ln(l/Z) 



(25) 



The normalized friction coefficient (Eq. (19)) is thus dominated by the second term and 
increases with increasing 

In the opposite case, r>l (T± ^> T\\), the evoluation of Eqs. (20) and (21) yields 



h 
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mo 



l7¥-l-21ni±Vi±ir 
2 2 



7T 



mo 



1 + 



i+ir 



(26) 



(27) 
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and 



h = ~ h * 2zmr> (28) 

Then h+h sin 2 ■§ ~ ii cos 2 d and the normalized friction coefficient decreases with increasing 
of d in this case. 

In Fig. 1 the normalized friction coefficient I\ + I2 sin 2 1? is plotted as a function of 
temperature anisotropy r for d = (solid line), d = n/Q (dotted line), $ = ir/3 (dashed 
line), $ = n/2 (dot-dashed line) and for fixed plasma density and average temperature 
(Z = 0.2). Fig. 1 shows an enhancement of the friction coefficient when the ion moves along 
the direction with low temperature. This effect can be easily explained in a binary collision 
picture. Let us consider the particular case of strongly anisotropic plasma T± >> Tj|. In this 
case the plasma electrons move mostly in the direction across to the anisotropy axis. For 
$ ~ 7r/2 the projectile ion moves along the plasma electrons thermal fluctuation direction 
and effective impact parameter for electron-ion collision is very small. Then the friction 
coefficient decreases. For ~ the projectile ion moves across to the direction of plasma 
electrons thermal fluctuation. Therefore, the impact parameter for electron-ion collisions 
increases which rises the friction coefficient. 

For arbitrary projectile velocities we evaluated Eq. (16) numerically. In Figs. 2 and 
3 the stopping power is plotted for strongly temperature anisotropic plasmas (r = 1CT 2 
and r = 10 2 in Figs. (2) and (3) respectively) with n = 10 8 cm -3 , T = O.leV and for 
four values of 1?; ■& — (dotted line), $ = n/Q (dashed line), = n/3 (long-dashed line), 
$ = 7r/2 (dot-dashed line). The solid lines are plotted for temperature isotropic plasma with 
T = T = O.leV. The general behaviour of the stopping power for two anisotropy parameters 
r is characterized by an increase by comparision with the isotropic case. At ~ 7r/2 and 
r = 10~ 2 (Fig. (2)) the ion moves in direction accross to the longitudinal electron motion 
with the lower temperature T± and the maximum of the stopping power is around V ~ fth±, 
whereas the maximum for an ion motion in longitudinal direction is at V ~ i> t h|| fth±- 
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IV STOPPING IN PLASMAS WITH WEAK MAGNETIC FIELD 



For the case when the magnetic field is weak, in the sense that the dimensionless parame- 
ter r) = oJ c /ujp is much less than unity, the functions G and F, Eqs. (6) and (7), which define 
the dielectric function, can be expanded about its field free values go{s/A)/A 2 , fo(s/A)/A 2 
Eqs. (14) and (15) 



G(s)+zF(s) = ±- 2 r 



90 (£) + %h {j 



+ V 2 ^^[9i(s)+ i f 1 (s)}, (29) 



where 



2 7 ft 2 \ 

9i(s)+ifi(s) = g(!- r ) Jt 3 dtl—Tsm 2 a-ljexp(istV2-A 2 t 2 

r— OO 

+ —— r [ t 4 dtexp(istV2 - A 2 t 2 ), 
6 J 
o 

s = u>/kvth.\\. Substituting this expression (29) into Eq. (8) leads to 



(30) 



S = S + V 2 S U (31) 

where So is the stopping power in plasma without magnetic field Eq. (16) and ifSi represents 
the change due to a weak magnetic field. After some simplifications it becomes 



rr Z 2 e 2 V }, J, (1 - /x 2 ) cos 2 9 



^2 2^^/ ^ 7. (32) 

X GXp I I 

I U 2 . I, 2A 2 J f 2 ( V cos@ \ , n 2 ( V cos 9 

v 11 ; jo (v^s-j +9 ° 

In the temperature isotropic plasma (r = 1) Eq. (32) coincides with the results by May and 
Cramer J?| after integration over ip. Note that the additional term Si does not depend on 
the cutoff parameter fc max - 

In the next subsections we evaluate Eq. (32) for small and large projectile velocities. 
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A Small projectile velocities 



When the projectile ion moves slowly (V < v th ) in plasma Eq. (32) leads to the simplified 
expression 



,S', = ^U/-^— P(&,t), (33) 
607T A 2 n V 2 v th y h 



with 



/ 1 _|_ T \ 3 / 2 r 

P(#,T)=(^f-) [j>,(r) + J> 2 (T) S m 2 tf], (34) 

p ^=w^\ Ut+25 - 3J j^ Mt)} ' (35) 

W = T^ibp {^^T=T Mt) - 28T ' 2 ~ 91t + 2 > ' (36) 

Here, the function Pq{t) is given by Eq. (12). In temperature isotropic plasma with r = 1 
we have P x (l) = P 2 (l) = 1. 

In Fig. 4 the normalized friction coefficient P("#, r) for the additional stopping power Si 
is plotted as a function of r for = (solid line), $ = n/6 (dotted line), = it/ 3 (dashed 
line), $ = tx/2 (dot-dashed line). The general behaviour of P(#, r) is similar to the friction 
coefficient of the plasma without magnetic field (see Fig. 1). Here, the correction P($, r) 
can be also negative at small r and which then corresponds to a slight decrease of the 
stopping power, Eq. (31). 



B High projectile velocities 



When the projectile ion moves with large velocity {V ^> i> t h), Eq. (32), yields 



Z 2 e 2 ul 
8V 2 



' 2 f 

^ 2d(l + cos 2 tf) - C 2 P(tf, r) 



cos 2 d + 



sin i? 
B(i?,r) + 1_ 



(37) 
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where 



/ t \ V2 

mr) = [ 2 o ■ ■ 2 J » (38) 

V t cos z v + sin 1/ / 

1 7x 2 exp(— x 2 /2)dx 1 7x 2 (7 — x 2 ) exp(— x 2 /2)dx 

1 = 3^/ /o 2 (x)+^ 2 (^) ' = 37^/ Zo 2 (*)+SoV) ' ( } 

For numbers Ci and C2 we get the accurate values C\ — 1 and C2 = (C2 — 10~ 12 ) 

respectively. Therefore from Eq. (37) we have finally 

Si = -^^(1 + cos 2 0). (40) 

This result is in accord with the results of Honda et al. || and May and Cramer 0, who, 
however, kept the terms 0(V~ 4 ) in their work as well. Although the function Si in Eq. (40) 
is proportional to the plasma density, the full correction term ^ 2 Si does not depend on the 
plasma density. 

In Figs. (5) and (6) we show the velocity dependence of the function Si for r = 10~ 2 
and r = 10 2 respectively. The different curves are $ = (solid line), = n/6 (dotted line), 
$ = 7r/3 (dashed line), d = it /2 (dot-dashed line). For small and medium projectile velocities 
the weak magnetic field decreases the total stopping power for small r and increases it in the 
high t limit. For high projectile velocities the magnetic field always reduces the stopping 
power independent of the temperature anisotropy, see Eq. (40). 



V STOPPING IN PLASMAS WITH STRONG MAGNETIC FIELD 

We now turn to the case when a projectile ion moves in a temperature anisotropic plasma 
with a strong magnetic field, which is on one hand, sufficiently weak to allow a classical 
description (hu c < k^Tj_ or h/ m%i < a c , and, on the other hand, comparatively strong so 
that the cyclotron frequency of the plasma electrons exceeds the plasma frequency u c ^> uo p . 
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This limits the values of the magnetic field itself and values of perpendicular temperature 
and plasma density. From these conditions we can obtain 



3 x KT 6 nJ /2 <B < 10 5 T ± , (41) 

where n is measured in cm -3 , T± is measured in eV, and B in kG. Conditions (41) are 
always true in the range of parameters no < 10 15 cm~ 3 , Bq < lOOkG, T± > 10~ 3 eV. Then the 
perpendicular motion of the electrons is completely quenched and the stopping power de- 
pends only on the longitudinal electron temperature Tji . The dependence on the transversal 
temperature will be only introduced by the cutoff parameter Eq. (9). 
In the limit of sufficiently strong magnetic field, Eq. (8) becomes 

2 ^ 2 /\ 3 „ h h cose/oOO 



Sinf = 2T2 k dk dp dip—— 2 (42) 

* X d\\{ { { [k 2 + 9o(s)] +fi(s) 

with s = (K/f t h||)(cos 0//x) and g , fo from Eqs. (14), which gives after integration over k 



Z 2 e 2 r } ( V cos 6 \ 
5'inf = - 2x2 dn dip cos6<2o ,£|| • (43) 



2n 2 X 2 m J J ' \v tH « 
Here the function Q is given by Eq. (17). For further simplification of Eq. (43) we introduce 
the new variable of integration x = cosB/u. After p integration in Eq. (43) we finally find 
the stopping power in the presence of a strong magnetic field as 



Z 2 e 2 „ ( V 

where 



Sm(Y,-&) = —- —Q[ ,#), (44) 



Q —,$) = sin 2 tf / , 9 v ""' / —— ^. (45) 

\v t h\\ / J Or 2 + i-2x cos??) 3 / 2 v ' 



In the previous works [11, 12, Q3[ only the case of ■& = the motion of the projectile 



along the magnetic field direction has been investigated. In this case the integral in Eq. (45) 
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diverges, while prefactor sin 2 1? tends to zero. Introducing the new variable of the integration 
in Eq. (45) y = (x — cost?)/ sin$ we obtain for vanishing angle $ 



Q 



v 



■d -> 



2Q 



V 



en 



(46) 



V t h\\ J \Vth\\ 

Thus expression (44) reproduces the known results for the stopping power on an ion 



which moves along the direction of the magnetic field [11, 12 , [13| . 

In the following paragraphs we will discuss its low and high velocity limits. 



A Small projectile velocities 



In the low velocity limit (V -C v t h\\) Eq. (45) becomes 



Q ( — ,0 



V t h\\ 



2V 

Vth\\ 



where 



2^(^11 



sin -#ln 



2v 



V siwd 



+ 1-2 sin 2 i? 



+ C 1 (en)sin 2 ^}, (47) 



Ci(£||) = / % [Qo(*,£||) - V^(£||)sl + [^Qo(x,Z\\). 



(48) 

1 

Here, the function ip is defined by Eq. (22). Since we deal with small ion beam-plasma 
coupling 2 < 1 we have, £y 3> 1 in Eqs. (47) and (48) and the function Ci(£) simplifies 



Ci(Cii) ^ v^Fln-ln^|| +0.6, 

1 T 



(49) 



where 7 = 0.5772 is Euler's constant. 

We note that the friction coefficient S-^/V from Eqs. (44) and (47) contains a logarith- 
mically large term which vanishes for $ — > 0. It will be shown in the next section that this 
behaviour is a characteristic feature of the stopping power at low velocities and the friction 
coefficient for arbitrary strength of the magnetic field. 
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B High projectile velocities 



In the case of high projectile velocities (V 3> v t h\\) the general expression (45) becomes 



Q^U^Isin^ 



«th|| 



V 2 



ln ( 2V J + C 2 (£n 



(50) 



where 



i 1 00 j 

1 



Ztc 



(51) 



which gives for £ii ^> 1 C*2 (C|| ) — l n £||- The stopping power for strong magnetic fields shows 
in the low and high velocity limits (Eqs. (47) and (50)) an enhancement for ions moving 
transversal to the magnetic field compared to the case of the longitudinal motion ($ = 0). 
This effect is in agreement with PIC simulation results |TJ]]. In contrast to the field-free 
case, at strong magnetic field and for $ = 0, V ^> v th \\ (Eqs. (44) and (50)) we have 
S in { ~ Z 2 e 2 u 2 /2V 2 independent of A; max . The cutoff /c max necessary at low ion velocities is, 
however, less well defined here than for the field-free case, where the cutoff (9) was deduced 
from the binary collision picture. Now, the electrons are forced to move parallel to Bo- Since 
we assumed the motion of the ion in this direction as well the ion and an electron just pass 
each other along a straight line. For symmetry reasons the total momentum transfer and 
the stopping power is zero. Purely binary interactions contribute nothing and the stopping 
of the ion is only due to the collective response of the plasma, that is, due to modes with 
long wavelengths k < 1/\d\\- This suggests taking k m3iX of the order of 1/\d\\, but further 
investigations are clearly needed here for a more precise description in this particular case. 

In Figs. (7) and (8), the stopping power S^f is plotted as a function of projectile velocity 
(in units of v th{{ ) for n = 10 6 cm" 3 , 7], = 10~ 4 eV , T ± = 10" 5 eV (Fig. (7)), T ± = O.leV 
(Fig. (8)), and for four different values of angle i? : i? = (solid line), $ = ir/6 (dotted 
line), $ = 7r/3 (dashed line) and = 7r/2 (dash-dotted line). The enhancement of Si n f(V, $) 
with respect to Sinf (V, 0) in the low and in high velocity limit by increasing of the angle $ is 
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documented in Fig. (9), for T\\ = l(T 4 eV, T± = O.leV, n = 10 6 cm 3 , $ = 7r/6 (solid line), 
i? = 7r/4 (dotted line), i? = 7r/3 (dashed line) and i? = 7r/2 (dash-dotted line). The physical 
origin of this angular behaviour in the low and high velocity limits is the enhancement of 
the effective impact parameter for an individual electron-ion collision with increasing For 
medium projectile velocities V ~ v t M\ the collective excitations in plasma become important 
and then stopping power is higher for small 

VI STOPPING AT ARBITRARY MAGNETIC FIELD AND IN 
LOW- VELOCITY LIMIT. ANOMALOUS FRICTION COEFFICIENT 

We now proceed with a projectile ion at low velocities and at arbitrary magnetic field. 
This regime is of particular importance for the electron cooling process 0, 0, 0]- In the 
presence of a magnetic field the friction coefficient here contains a term which diverges like 
\n(vth\\/V) in addition to the usual (see e.g. Sec. Ill) constant one. 

For this consideration it is convenient to use the Bessel function representation of the 
dielectric function which has been given e.g. by Ichimaru [fnj] , see Appendix A Eq. (A7), 
and to write the real and imaginary parts of Eq. (A7) separately 
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The notations in Eqs. (52) and (53) are explained in Appendix A. 

For the friction coefficient we have to consider S, given by Eq. (8) in the low-velocity 
limit and thus the functions G and F given by Eqs. (52) and (53), when oo = kV. Now we 
have to write the Taylor expansion of Eqs. (52) and (53) for small oo = kV. However, the 
first term of Eq. (53) exhibits a singular behaviour in the limit of u — kV — > where the k\\ 
integration diverges logarithmically for small fcy. We must therefore keep u = kV finite in 
that integration to avoid such a divergence. This anomalous contribution which arises from 
the first term of Eq. (53) in low-velocity limit is 
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"an — 



2 X 1 / 2 Z 2 e 2 V 
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A (z) exp 



v 2 cos 2 e 
2 < h |i ^ 
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where A (z) = exp(—z)I (z) and E 2 (k, ji) = G(oo = 0) is 



(54) 



«*">- 1 + 2 £(H£-M4»(i&)- 



(55) 



Here z = (k 2 r/r] 2 )(l — /i 2 ), /i = cos a = and is the angle between k and V. After \i 

and ip integration, see Appendix B, Eq. (54) reads 
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The function and thus S^n (56) vanishes in the limit B — > (or 77 — > 0) like 



(56) 
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The anomalous term Eqs. (56) and (57) therefore represents a new effect arising from 
the presence of the magnetic field, which is not restricted to anisotropic plasmas. 

For temperature isotropic plasma (r = 1) and for a sufficiently weak magnetic field rj < £y 
(or uj c < fcmaxfthii), Eq. (57) takes the form 

jF(r, 77) ~ exp 

where K and K\ are the modified Bessel functions of the second kind. In the case of very 
strong magnetic field 77 > i\\\fr (or uo c > k milx \ D ±) the function jF(r, rj) reads 

^(r, r,) ~ = ln(l + ) - • (60) 

The physical origin of such an anomalous friction coefficient may be traced to the spiral 
motion of the electrons along the magnetic field lines. These electrons naturally tend to 
couple strongly with long- wavelength fluctuations (i.e., small k\\) along the magnetic field. 
In addition, when such fluctuations are characterized by slow variation in time (i.e., small 
uj = kV), the contact time or the rate of energy exchange between the electrons and the 
fluctuations will be further enhanced. In a plasma such low-frequency fluctuations are pro- 
vided by the slow projectile ion. The above coupling can therefore be an efficient mechanism 
of energy exchange between the electrons and the projectile ion. In the limit of V — > 0, the 
frequency uj = kV — > tends to zero as well. The contact time thus becomes infinite and 
the friction coefficient diverges. 

The anomalous friction coefficient (see Eq. (56)) vanishes, however, when the ion moves 
along the magnetic field (i? = 0). Then the friction coefficient is solely given by the second 
term of Eq. (53). The contribution of this term to the stopping power leads to the usual 
friction law in plasma and reads for arbitrary angles $ 




X 2 m v tH J J fj, [k 2 + E 2 {k,^)} 2 



7T 



(61) 
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with 



Ei(k,jj) = ^n{z) exp 



n=l 



oc 



( 



2k 2 (i 2 




) 




(62) 



and E 2 (k,fi) as denned by Eq. (55). 

In Figs. (10) and (11) we compare the anomalous term S an with the low velocity stopping 
without magnetic field So see Eq. (19), where S an /S is plotted as a function of uo c /uo p for 
= tt/6 (solid line), d = tt/3 (dotted line), d = ir/2 (dashed line), Z = 0.1, V/v th = 0.2, 
and for two values of the anisotropy parameter r: r = 0.1 (Fig. (10)), r = 10 (Fig. (11)). 

We conclude that the anomalous term S an gives espessially for strong magnetic fields 
{uj c > u p ) and for strongly temperature anisotropic plasma (T± ^> Ty) an important contri- 
bution to the stopping. It should be noted that the observed enhancement of stopping due 
to S^n for Tj_ ^> T\\ can be potentially interesting for future electron cooling experiments. 



The purpose of this work was to analyze the stopping power of an ion in temperature 
anisotropic magnetized classical plasma. A general expression obtained for stopping power 
was analyzed in four particular cases: in a plasma without magnetic field; in a plasma with 
weak and very strong magnetic fields; and in a plasma with arbitrary magnetic field and for 
low-velocity projectile. 

From the results obtained in Sees. III-V, we found that the stopping power essentially 
depends on the plasma temperature anisotropy. In field-free case and for small ion velocities 
the effect of the anisotropy results in an enhancement of the stopping power when the ion 
moves in the direction with low temperature. 

For small projectile velocities a weak magnetic field slightly decreases the field-free stop- 
ping power for small r, in the opposite case (large r) the field-free stopping power slightly 
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increases. In the high-velocity limit correction to the field-free stopping power for weak 
magnetic fields is always negative and the stopping power is reduced by the magnetic field. 

In the case of strong magnetic fields we demonstrated an enhancement of the stopping 
power with increasing of d for low and high-velocity regions compared to the case of an ion 
which moves along B . 

In low-velocity limit but for arbitrary magnetic field, we find an enhanced stopping 
power compared to the field-free value mainly due to the strong coupling between the spiral 
motion of the electrons and the long-wavelength, low-frequency fluctuations excited by the 
projectile ion. This anomalous stopping power increases with the angle $ (the angle between 
ion velocity V and magnetic field B ) and depends strongly on the temperature anisotropy 
t = T±/T\\, as seen in Figs. (10) and (11). Although the nature of the anomalous stopping 
power is conditioned by the external magnetic field the temperature anisotropy of the plasma 
can intensify this effect when T± 3> Ty (see Fig. (11)). 

This emphasizes the importance of the special role of fluctuations with small k\\ and small 
uj (small projectile velocity V) and as another significant contribution to the energy exchange 
processes arising from the collective modes of plasma. Potentially, the electron plasma waves 
and the ion acoustic waves in a magnetized plasma might provide a significant energy- 
exchange mechanism between projectile ion and plasma particles. This fact then makes 
it necessary to consider the influence of plasma collective modes to anomalous stopping 
process. This problem will be treated in a subsequent work. 
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A 



Here we describe the evaluation of the dielectric function in the temperature anisotropic 
case where the velocity distribution of the unperturbated distribution function was given 
by Eq. (4). We next introduce the Fourier transformations of /i(r, v, t) with respect to 
variables r and t, /i(k, 00, v). Because of the cylindrical symmetry (around the magnetic 
field direction b = B0/-B0 — z) of the problem, we choose 



v = v± cos ax + v± sin ay + v\\z. 
Then the Vlasov Eq. (2) for the distribution function becomes 



(1) 



d i 

— fi(k,u, v) H (kv -u- iO)/i(k, u, v) 

OCT Ul r 



— — (k^ 

mo;,, \ <7v 



(2) 



where </?(k, u;) is the Fourier transformation of tp(r,t). The positive infinitesimal +i0 in Eq. 
(A2) serves to assure the adiabatic turning on of the disturbance and thereby to guarantee 
the causality of the response. The solution of the Eq. (A2) has the form 



/i(k,o;,v) = ip{k,u) da 2 [k—) exp — / da x [-u - iO + (kv) CT: 



Combining Eq. (A3) with the Poisson equation (3) we find for the dielectric function 
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(4) 



where k x = k± cos ky — A^j_ siny?. After integration by the variables a± : a 2 and a, and 



using the expression 2C 



exp(— iz sin 0) = V] J n (z) exp(—in9) , 



(5) 
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where J n is the Bessel function of the nth order, we obtain the expression 



g7r 2 e 2 +00 - Y ( 

e(k,u;) = l-— £ jv ± dv ± J dvA 

n— 00 _ QO \ 



' nujcdfp {k d fo\ Jn( k ± V ±/ UJ c) 
fj_ " y nu; c + — c<j — iO 



Substituting Eq. (4) for the unperturbed distribution function f into Eq. (A6) we 
finally results in 



w 



u — nu r 



An(/3) , (7) 



lfclll^thll / 

where (3 = kj_v^ h± /uj^ = k\a 2 c) A n (z) = exp(—z)I n (z), I n {z) is the modified Bessel function 
of the nth order, and W(z) is the plasma dispersion function | jT%|| . 

To show the identity of the two forms (Eqs. (6) and (A7)) of the dielectric function we 
will use the expansion in modified Bessel functions p0| 



exp(zcos9)= ^2 I n (z) exp(in9) . (8) 

n=— 00 

This allows to rewrite exp[— X(t)] with X(t) from Eq. (7) as 

exp[-X(t)} = exp(-i 2 cos 2 a) £ A n (/3) exp ffH^E] . (9) 

n=-oo \ kV th\\ / 

Substituting Eq. (A9) into expression (6) and integration over the variable t leads to 
Eq. (A7). 
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We now give a more detail derivation of the anomalous term S an (Eq. (56)). We start 
with the expression (see Eq. (54)) 

C^.A)=/^*(M,^xp(-^^), (1) 
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where <f>(fj,,(p) = cos 6, A = V/v th \\, 
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[k* + E 2 (k,fj.)Y 
For A — > a leading-term approximation of (Bl) leads to 



(2) 



Q(k,<p,\) ~ $(0,fc,p)ln- 



V2 



A|0(O,^)|V7 

where 7 is the Euler's constant, |0(O,y?)| = sm-&\ cos<^|, 
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Using the relation |17|, ^ 
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the function E2(k, 0) we finally takes the form 



(6) 



£ 2 ( M ) = 1 + (l - 1) K Q {k 2 r/rf). 



(7) 



Substituting Eqs. (B3), (B4) and (B7) into Eq. (54) and integration over cp we finally 
come to expression (56). 
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